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ABSTRACT 

Kepler- 16 is an eccentric low-mass eclipsing binary with a circumbinary transiting planet. Here we inves- 
tigate the angular momentum of the primary star, based on Kepler photometry and Keck spectroscopy. The 
primary star's rotation period is 35 . 1 ± 1 .0 days, and its projected obliquity with respect to the stellar binary 
orbit is 1.6 ±2.4 degrees. Therefore the three largest sources of angular momentum — the stellar orbit, the 
planetary orbit, and the primary's rotation — are all closely aligned. This finding supports a formation scenario 
involving accretion from a single disk. Alternatively, tides may have realigned the stars despite their relatively 
wide separation (0.2 AU), a hypothesis that is supported by the agreement between the measured rotation pe- 
riod and the "pseudosynchronous" period of tidal evolution theory. The rotation period, chromospheric activity 
level, and fractional light variations suggest a main-sequence age of 2-4 Gyr. Evolutionary models of low-mass 
stars can match the observed masses and radii of the primary and secondary stars to within about 3%. 

Subject headings: stars: binaries, rotation, activity, late-type, low-mass, formation, individual (Kepler-16 A, 
KIC 12644769) — planets and satellites: formation 



I. INTRODUCTION 

Kepler-16 (AB)-b is a planet with two parent stars (Doyle et 
al. 201 1). The stars (0.7 and 0.2 M©) are in 41-day eccentric 
orbit, and the planet (0.3 Mjup) circles both of them every 229 
days. Viewed from the Solar system, the stars eclipse each 
other and the planet transits both of them, providing definitive 
evidence that circumbinary planets exist and permitting pre- 
cise determinations of the system's parameters. For example 
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the planet's radius is known to within 0.3%, better than that 
of any other known exoplanet. The stars are themselves of 
interest as a rare example of low-mass dwarfs with precisely 
known dimensions. 

Such a unique system should be studied in every possible 
way, for exploratory purposes as well as the specific purpose 
of understanding its formation and evolution. How old are the 
stars? Did the planet form together with the stars, or was it 
captured from another system? Has there been tidal evolution 
or other effects that have modified the system's architecture? 
Here we present an investigation of the angular momentum of 
the primary star, bearing on these questions. 

It has already been established that the planes of the cir- 
cumbinary orbit and the stellar orbit are aligned to within 
0.°5 (Doyle et al. 2011). This suggests all three bodies inher- 
ited their angular momentum from a single disk, as opposed 
to dynamical scenarios that are often invoked for triple sys- 
tems such as close encounters (Mikkola 1984, Bailyn 1989, 
Ivanova 2008) or dynamical decay (Sterzik & Tokovinin 
2002). One must remember, though, that the planet was dis- 
covered with transit photometry, a technique that is severely 
biased toward finding coplanar orbits. This raises the ques- 
tion of whether the orbital coplanarity of Kepler-16 is at all 
representative of circumbinary planets, and motivates mea- 
surements of the alignment between the orbital axes and the 
stellar spin axes, for which there was no selection bias. 

This Letter is organized as follows. Section 2 presents a 
photometric determination of the rotation period. Section 3 
presents a spectroscopic determination of the sky-projected 
stellar obliquity (the angle between the rotation axis of the 
primary star, and the stellar orbital axis), based on observa- 
tions of the Rossiter-McLaughlin effect. Section 4 discusses 
the implications of these results for our understanding of the 
primary star and of the system's history. 
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2. THE ROTATION PERIOD 

We measured the rotation period of the primary star with 

data from the Kepler spacecraft, a 0.95m space telescope 
that monitors the optical brightness of about 150,000 stars 
in a quest to detect transits of potentially habitable Earth- 
sized planets (Borucki et al. 2010). Overviews of the mis- 
sion design, the instrument performance, and the data pro- 
cessing pipeline were given by Koch et al. (2010), CaldweU 
et al. (2010) and Jenkins et al. (2010). 

Kepler-16 was observed with 29.4 min sampling for a 
nearly continuous 600-day interval, from 2009 May 02 to 
2010 December 22 (quarters 1-7). The duty cycle was 94%, 
with 17 short gaps due to technical problems as well as sched- 
uled interruptions in observing. After each interruption a 
jump was observed in the relative flux. We placed all the data 
onto a common flux scale under the assumption that the flux 
variations during the interruptions were smooth enough to be 
described by a quadratic function of time. Specifically, we 
multiplied the data from each of the 18 disjoint intervals by a 
constant, and determined the optimal values of the constants 
by fitting quadratic functions to the data within one day of an 
interruption. 

The resulting time series exhibited a secular 3% decrease 
in relative flux, which could be an instrumental effect or a 
true decrease in steUar brightness. Since this ttend is irrele- 
vant to the rotation period determination, we applied a 70-day 
median filter prior to plotting the time series in the top panel 
of Figure 1 . The time series exhibits quasiperiodic variations 
of order 0.5%. As usual for late-type stars, we attribute these 
variations to dark spots and bright plages being carried around 
by stellar rotation. 

The bottom panel of Figure 1 is a Lomb-Scargle peri- 
odogram, showing a prominent peak at 35.1 days along with 
smaller peaks at the first two harmonics. We identified this 
peak with the stellar rotation period. We estimated the uncer- 
tainty in the period by dividing the data chronologically into 4 
equal segments, analyzing each piece separately, and finding 
the standard deviation in the mean of the periodogram peaks. 
Based on this analysis we find Prot = 35 . 1 ± 1 .0 days. 

3. THE ROSSITER-MCLAUGHLIN EFFECT 

We measured the sky-projected obliquity and rotation rate 
of the primary star by conducting spectroscopic observations 
of a primary eclipse and analyzing the Rossiter-McLaughlin 
(RM) effect. The RM effect is the anomalous Doppler shift 
that is observed during eclipses as a consequence of the se- 
lective blockage of the rotating stellar photosphere (Rossiter 
1924, McLaughHn 1924). 

We used the Keck 1 10m telescope and HIRES spectrograph 
to gather 14 spectra on 2011 May 28/29, starting 40 min- 
utes before ingress and extending for 5 hours until morning 
twilight, thereby covering about three-quarters of the eclipse. 
Another 3 spectra were obtained the following night to track 
the out-of-eclipse velocity variation. The typical exposure 
time was 19 minutes. The I2 absorption cell was used to 
establish the wavelength scale and instrumental profile. A 
single exposure without I2 was also obtained to serve as a 
template spectrum. The relative radial velocities (RVs) were 
determined with a descendant of the algorithm of Butler et 
al. (1996). They are given in Table 1 and plotted in Figure 2. 

The "red-then-blue" pattern of the anomalous Doppler shift 
is characteristic of a prograde orbit with good aUgnment be- 
tween the primary star's rotational and orbital angular mo- 



menta. In the first half of the eclipse, the secondary cov- 
ers the approaching (blue) half of the primary, causing the 
net starUght to be redshifted; then, the secondary moves over 
the receding (redshifted) half of the primary, producing an 
anomalous blueshift. For quantitative modeling we used the 
technique of Albrecht et al. (2007), in which a pixellated stel- 
lar disk is constructed, and a theoretical spectral line profile is 
computed for each pixel based on the local intensity and ve- 
locity of the photosphere. The integrated spectrum is obtained 
by summing over the unecUpsed pixels, and then the RV is 
calculated by cross-correlation with the integrated spectrum. 

To compute the relative intensities of the pixels, we as- 
sumed a linear limb-darkening law. The pixel velocities in- 
cluded the effects of uniform rotation, macroturbulence, and 
the convective blueshift. The model for macroturbulence was 
taken from Gray (2005), assuming equal radial and tangen- 
tial velocity perturbations with a standard deviation C,ky = 
2.5 km s~' . The model for the convective blueshift was taken 
from Shporer & Brown (201 1), in which the velocity of each 
pixel is shifted by Vcb = 0.2 km s"' away from the center of the 
star. The stellar radii and orbital inclination were held fixed 
at the values determined by Doyle et al. (201 1). We neglected 
any light from the secondary, as the light ratio is constrained 
to be <1.3% (Doyle et al. 2011). 

The 6 adjustable parameters were the projected rotation rate 
vsinj, the projected obliquity (5, the limb darkening coeffi- 
cient u, the central ecUpse time Tc, the zero point of the rela- 
tive RV scale 7, and a parameter describing the out-of-ecUpse 
RV variation. For the latter, a simple linear slope would have 
sufficed, but for convenience we used M'g, the mass of the sec- 
ondary when all other orbital parameters (period, inclination, 
eccentricity, argument of pericenter) are held fixed at the best- 
fitting values reported by Doyle et al. (201 1). We expect small 
differences between M'^ and the true secondary mass Mb due 
to the small uncertainties in the other orbital parameters as 
well as spurious radial accelerations caused by starspots or 
Ught from the secondary. 

To determine the allowed parameter ranges we used 
a Markov Chain Monte Carlo algorithm, with the the 
Metropolis-Hastings algorithm and Gibbs sampler. The like- 
lihood was taken to be exp(-x^/2), where is the usual sum 
of the standardized residuals between the observed and calcu- 
lated RVs. Uniform priors were adopted for vsin;, (3, 7, and 
Mj. A Gaussian prior was used for Tc, with a central value 
equal to the predicted ephemeris time and a standard devia- 
tion of 1 min (the typical level of eclipse timing variations). 
A Gaussian prior was used for the limb darkening coefficient, 
M = 0.8 ±0.1, based on the tables of Claret (2000). 

The formal "lu" uncertainty interval was taken to be the 
range between the 15.8% and 84.2% levels of the cumulative 
distribution of the marginalized posterior for each parameter 
In addition, we checked on the sensitivity of the results to the 
assumed values of the macroturbulent velocity (^rt and con- 
vective blueshift Vcb, by perturbing each of those quantities 
by ±50% and refitting. (It would be better to allow these pa- 
rameters to vary during the fit, but in practice this was too 
computationally demanding.) Changing the macroturbulence 
had no appreciable effect. Changing the convective blueshift 
caused shifts of 18 m s~' in vsinj and 2° in (3, comparable to 
the formal Icr intervals. Consequently, we enlarged the un- 
certainty intervals for those two parameters by adding those 
shifts in quadrature with the formal errors. Table 2 gives the 
results for all the parameters. 
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Fig. 1. — Kepler photometry. Top. — Time series of relative fiux. Eclipse data have been removed. Vertical dotted lines are plotted every 35.1 days, the 
estimated rotation period. Different colors indicate the 4 segments for which periodograms were computed separately in order to gauge the uncertainty in the 
rotation period. Bottom. — Lomb-Scargle periodograms of the entire time series (thick black line) and each of the 4 segments (thinner colored lines). 



4. DISCUSSION 
4.1. Spin-orbit alignment 

The angle (5 between the sky projections of the primary's 
rotational angular momentum and the orbital angular momen- 
tum of the binary was found to be 1 .6 ± 2.4 degrees|]3 The star 
is apparently aligned with the orbit to within a few degrees. 

Ordinarily a warning must be issued here: the stellar ro- 
tation axis might be inclined along the line of sight, and the 
sky-projected obliquity [3 might not be representative of the 
true three-dimensional obliquity -ip. In this case, though, an 
upper bound on the true obliquity is enforced by the com- 
bination vsin/ = 0.920 ±0.025 km s"' from the RM anal- 
ysis, P,ot = 35.1 ± 1.0 d from the Kepler data, and R^, = 
0.6489 ± 0.0013 Rq from the photometric-dynamical model 
of Doyle et al. (2011). Assuming v = lirRi,/ P^ot, and adopt- 
ing a isotropic prior for / (uniform in cos/), we find sin / = 
0.994!|;{j°f and / = 90° ± 9°. With 95.4% confidence, the true 
obhquity is < 18?3. 

This system's angular momentum has 5 contributions — 
the stellar orbit, the planetary orbit, the primary rotation, the 
secondary rotation, and the planetary rotation — with magni- 
tudes in the approximate ratios 10000 : 40 : 1 : 0.1 : 0.001. 
Kepler- 16 is an orderly system, with good alignment between 
the three largest portions of the angular momentum. To our 
knowledge, Kepler- 16 is the longest-period stellar binary for 
which a stellar obliquity has been measured (see Table 1 of 
Albrecht etal. 2011). 



The angle /3 is defined using the coordinate system of Hosokawa (1953). 
Note that Ohta et al. (2005) and others have used a different symbol A to 
represent this angle, and a coordinate system such that A = — /3. 



4.2. Activity, rotation, and age 

The Keck spectra show Ca H&K chromospheric emission 
with log/?Jjj^ = -4.68 ±0.10, stronger than the Sun's value of 
-4.91. The spectra obtained previously by Doyle et al. (201 1) 
also exhibit Ca H&K emission, with RV variations tracking 
those of the primary, proving that the emission originates on 
the primary rather than the secondary or the planet. We may 
therefore use empirical relations between the rotation, chro- 
mospheric activity, and age of main-sequence dwarfs to esti- 
mate the age of Kepler- 16 A and check whether there is any- 
thing unusual about its properties. 

Building on work by Barnes (2007) and others, Mamajek & 
Hillenbrand (2008) provided up-to-date activity/rotation/age 
relations forK2-F7 dwarfs (?aO.8-1.3M0), which we extrapo- 
lated to interpret the 0.7 Mq primary of Kepler-16. The activ- 
ity/age relation in their Eqn. (3) gives an age of 1 .8 ± 1 .2 Gyr 
The rotation/age relation implicit in their Eqns. (12-14) give 
an age of 3.7 ±0.8 Gyr A similar age is obtained from the 
relations of Barnes (2010) and Barnes & Kim (2010). 

In addition, the rotation and activity are expected to be di- 
rectly linked, with a particularly strong coiTelation betwen the 
Rossby number Ro (the ratio of rotation period to convective 
turnover timescale) and logTJjjj^. We used Eqn. (4) of Noyes 
et al. (1984) to estimate the convective turnover timescale, ob- 
taining Tc = 23.0 d and Ro = 1.52. The relation between Ro 
and chromospheric activity shown in Fig. 7 of Mamajek & 
Hillenbrand (2008) predicts log/?(jj^ = -4.70, in good agree- 
ment with the measured value. The light variations of 0.5- 
1 % are also typical for a star with this level of chromospheric 
emission (see, e.g.. Hall et al. 2009). 

All together, the rotation period, chromospheric emission 
level, and fractional light variations paint a picture of an ordi- 
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Fig. 2. — Keck radial velocities. Top. — Appai'ent radial velocity (solid points) and the best-fitting model (gray curve). Bottom. — After subtracting the 
best-fitting orbital model, thereby isolating the RM anomaly. Each night's data is shown separately, along with the residuals. 



nary 0.7 Mq dwarf star with an age of 2-4 Gyr 

4.3. Comparison to evolutionary models 

The Keck spectrum provides new estimates of the primary's 
photospheric parameters. Analysis with Spectroscopy Made 
Easy, a software package written by Valenti & Piskunov 
(1996), gives T,ff = 4337 ± 80 K and [Fe/H] = -0.04 ± 0.08 
for a fixed value of \ogg = 4.6527 (the value determined by 
Doyle et al. 2011, which has a negligible uncertainty for 
this purpose). These agree with the previously reported re- 
sults Teff = 4450 ± 150 K and [m/H] = -0.3 ±0.2 (Doyle et 
al. 2011). 

Figure 3 shows a comparison between the observed masses 
and radii of Kepler- 16 A and B and the theoretical evolution- 
ary models of Baraffe et al. (1998). Models are shown for 1 
and 5 Gyr for a mixing-length parameter of 1 .0. The 1 Gyr 
model gives the best match to the primary mass and radius. 



although it is also possible to obtain a good fit at 5 Gyr by in- 
creasing the mixing-length parameter toward the Solar value 
of 1.9 (as might be expected for the primary; see, e.g., De- 
mory et al. 2009). The models also match the measured ef- 
fective temperature of the primary. Interestingly the calcu- 
lated primary radius is within a few percent of the observed 
radius. This is in contrast to the other stars of similar mass 
shown in Figure 3, for which the observed radii are 10-15% 
larger than the model predictions. Such discrepancies have 
been attributed to high activity and rapid rotation in the stars 
that have been studied closely (see, e.g., Chabrier et al. 2007). 
The secondary star's mass and radius are also within a few 
percent of the calculated values. We leave a more detailed 
comparison with models for future work, which should take 
into account the uncertainty in the metallicity as well as the 
a-element abundances (which have not yet been measured). 
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Fig. 3. — Theoretical and observed masses and radii of low-mass stars. The model curves are from Baraffe et al. (1998), for a metallicity of -0.04 (obtained 
by linear interpolation between 0.00 and —0.50) and a mixing-length parameter qml = 10. In addition to Kepler-16 (Doyle et al. 2011), data are shown for other 
low-mass stars for which the masses and radii have been determined to within 3% according to the rigorous criteria used by ToiTes et al. (2010). Those systems 
are: KOI-126 (Carter et al. 201 1), YY Gem (Torres & Ribas 2002), CU Cnc (Ribas 2003), GU Boo (Lopez-Morales & Ribas 2005), NGC 2204 S892 (Rozyczka 
et al. 2009), and IM Vir (Morales et al. 2009). 



4.4. Formation scenarios 

Binary stars are thought to form by fragmentation of col- 
lapsing prestellar cores, with close binaries such as Kepler-16 
possibly resulting from fragmentation during a late, isother- 
mal phase of the collapse (see, e.g., Goodwin et al. 2007). 
Stable hierarchial triples can also be produced by fragmenta- 
tion (Boss 1991, Bate 2009), raising the question of whether 
the three bodies of Kepler-16 formed in this manner 

In this scenario one would expect a tendency for the frag- 
ments to be aligned, but there is no obvious reason why the 
alignment would be as close as is observed for Kepler-16. 
Among the triple systems that formed in a large-scale sim- 
ulation by Bate (2009), the orbital planes were typically mis- 
aligned by «60°. Furthermore there are some close binaries 
with much larger obliquities than Kepler-16 A, such as DI 
Herculis (Albrecht et al. 2009). 

This suggests an additional chapter is needed in the story, 
beyond fragmentation. Perhaps after their formation the 
stars continued to accrete substantially from a circumbinary 
disk, which would have decreased their orbital separation and 
aligned their spin axes (Bate et al. 2002). Or the primary could 
have formed with a massive circumstellar disk, which then 
fragmented to form the secondary. In either of these scenar- 
ios, the planet could have formed by core accretion near the 
inner edge of the circumbinary disk (see, e.g., Pierens & Nel- 



son 2008, Mai-zari et al. 2008). 

4.5. Evidence for tidal evolution 

Another process that could have reduced the stellar obliq- 
uity is tidal evolution. In the fullness of time, tides synchro- 
nize and align the spins of a binary system, and circularize 
their orbit (Zahn 1977, Hut 1981). Conventional wisdom 
would say that tides are irrelevant for Kepler-16 due to the 
relatively long orbital period (41 days) and wide separation 
{a = 0.22 AU, a/Ri, = 74), and indeed the orbital eccentric- 
ity of 0.16 shows that circularization is incomplete. On the 
other hand, synchronization and alignment should be faster 
than circularization because the rotational angular momentum 
is smaller than the orbital angular momentum. More broadly, 
tidal evolution timescales are poorly known, especially for 
spin evolution (see, e.g., Mazeh 2008). 

An order-of-magnitude assessment of the importance of 
tidal evolution begins with the observation that most late-type 
binaries with periods <10 days have circular orbits (Mazeh 
2008). Therefore the circularization timescale is ^5 Gyr 
for a/Ri, ^ 20. If this timescale varies as (a/R^,)^ (Zahn 
1977), then for Kepler-16 it is 5 (74/20)** - 2 x 10^ Gyr, 
i.e., consistent with the observation that the orbit is still ec- 
centric. We further suppose that the spin evolution timescale 
is smaller by a factor of ^lO'* (the ratio of orbital to rota- 
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tional angular momenta) and varies as (a/R^,)^ (Zahn 1977). 
Then the spin evolution timescale for Kepler- 16 would be 
5 X 10"^ (74/20)'' ~ 1 Gyr. This is corroborated by the more 
detailed tidal model of Terquem et al. (1998): using their 
Eqns. (41-42), the timescales for circularization and synchro- 
nization are 10"* Gyr and 2 Gyr, respectively. These calcula- 
tions are subject to the well-known uncertainties in tidal evo- 
lution timescales, but they do suggest that the low obliquity of 
Kepler-16 A is at least partly a consequence of tidal evolution. 

In this context, the rotation period of the primary is intrigu- 
ing. According to the tidal theory of Hut (1981), before the 
orbit circularizes the stellar spins evolve into a pseudosyn- 
chronous state, with spin periods shorter than the orbital pe- 
riod due to the enhanced tidal dissipation at pericenter pas- 
sages. For Kepler-16 A, using Eqn. (42) of Hut (1981) with 
e = 0.15944 (Doyle et al. 2011), the predicted pseudosyn- 
chronous period is 35.62 days, in agreement with the mea- 
sured period of 35.1 ± 1.0 days. As noted in Section 4.1, 
though, the measured rotation period is not unusual even 
for an isolated star, and the agreement with the pseudosyn- 



chronous value may be a coincidence. A more detailed study 
of the tidal evolution of this unique triple system is warranted, 
as are further spin-orbit studies of relatively wide binaries. 
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TABLE 1 

Relative Radial Velocity Measurements of Kepler-16 A 



BJDuTC 


RV [ms-'] 


Unc. [m s 


2455704.90651 


405.81 


1.88 


2455 /U4.y2Uo7 


385.56 


1.88 


2455704.93459 


360.90 


1.91 


2455704. 94s35 


353.70 


1.88 


2455 /04. 96150 


346.34 


2.02 


2455704.97574 


340.26 


1.72 


2455704.99387 


321.48 


1.88 


2455705.00766 


302.26 


1.94 


2455705.02716 


262.90 


1.88 


2455705.04061 


222.41 


1.94 


2455705.05985 


167.53 


1.87 


2455705.07461 


127.95 


2.39 


2455705.08883 


86.27 


1.88 


2455705.10192 


56.34 


2.05 


2455705.90834 


-1129.19 


2.78 


2455705.96133 


-1207.97 


2.86 


2455706.11341 


-1447.21 


2.64 



Table 2: Summary of Results 



Parameter Value 

Projected rotation speed, v sin ; [km s"' ] 0.920 ± 0.025 

Projected spin-orbit angle, /3 [deg] 1.6 ± 2.4 

Secondary mass parameter, [Mq] 0.20133 ± 0.00026 

Velocity zero point, f [m s"'] 427.6 ±1.3 

Limb darkening coefficient, u 0.717 ± 0.079 
Central ecUpse time, Tc [BJDurc] 2,455,705.05388 ± 0.00069 

Stellar rotation period, P,ot [days] 35.1 it 1.0 

Stellar inclination angle, i [deg] 90 ± 9 

Stellar obUquity, 4> [deg] <18.3 (95.4% conf.) 

Effective temperature. Tea [K] 4337 ± 80 

Iron abundance, [Fe/H] -0.04 ± 0.08 

Surface gravity, log(g [cm s^^)] 4.6527 ± 0.0017 

Main-sequence age [Gyr] 3.0 ± 1.0 



Note. — AH the stellar parameters refer to the primary star, Kepler-16 A, except for the secondary mass parameter M'^ (defined in the text). The result for 
M'g differs by 0.00120 ± 0.00071 Mq from the secondary mass Me = 0.20255 ± 0.00066 Mq determined by Doyle et al. (201 1). 



